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I GUT scale threshold corrections in minimal SU(5) supergravity grand 

^ [ unification are discussed. It is shown that predictions may be made despite 

I uncertainties associated with the high energy scale. A bound relating the 



strong coupling constant to the mass scales associated with proton decay and 



H ' supersymmetry is derived, and a sensitive probe of the underlying theory is 

outlined. In particular, low energy measurements can in principle determine 

the presence of Planck scale (1/Mpi) terms. 
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Over the past two decades, much attention has been given to the possibihty of unifying 
the three gauge groups of the Standard Model into one group. These GUT models are theo- 
retically more appealing than the Standard Model for various reasons The 1990 precision 
LEP data strongly indicated that supersymmetry is needed to achieve grand unification , 
and spurred on many new analyses of different aspects of SUSY GUTs. 

With the increasing precision of data, further predictions will be affected by supersymme- 
try threshold effects and the details of the model at high energy. In this note we examine the 
effects of the non-degeneracy of the super-heavy GUT spectrum, and the possible existence 
of non-renormalizable operators from Planck scale physics. A number of treatments of these 
issues have been given recently for the SU(5) model It was found that many of the 

predictions for low energy observables were blurred by the high scale effects. For example, it 
was argued that the SUSY scale cannot be determined by a more accurate measurement of 
as 0, and the rate of proton decay cannot be predicted from low energy data if additional 
Planck scale terms are present [^]. 

Nevertheless, we will show here that there are still predictions to be made in this model: 

a. The effect of Planck scale non-renormalizable terms becomes smaller as the value of 
0:3 is varied to lower values, so the lower limit on 0:3 is not lost. This is especially interesting 
since there is currently a disparity in the values of 0:3 between the measurements made 
at weak scale energies, and those made at lower energies [0. Resolving this disparity and 
refining the measured value of 0:3 will provide an important test of the model. 

h. Since Planck scale physics smears the correlation between 03 and the mass scale 
which governs proton decay, proton decay will be a sensitive probe of Planck scale physics. 
In particular it will be seen that by purely low-energy measurements one can determine 
experimentally the degree to which the dominant Planck scale term is present. Thus models 
of this type allow for the first time a test for the existence of Planck scale physics and 
whether Planck scale physics impinges on low-energy (electroweak scale) physics. 

With respect to the second item, progress has been made recently in deriving models 
similar to the ones considered here from string theory so that the gravitational smearing 
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may be calculable in principle. Proton decay would then be a sensitive test of string theory. 
Our model of GUT physics is defined by superpotential p!0| , p!l|] 



W = Ai[itr(S=^) + iMtr(S2)] 
+A2Hx(S^y + 3M' 5^)H^ 

+euvwxY H^M^^/iM^^ + H^M^^/sMy, (1) 

where S, H, and H form a 24, 5, and 5 of SU(5) respectively, M and M are 5 and 10 matter 
superfields, fi and /2 are Yukawa coupling constant matrices in the generation space, and the 
mass parameters M and M' are set equal to account for a light Higgs doublet. (This is a well 
known fine tuning problem with this model. We consider elsewhere alternate models such as 
those in Ref. [|1^] which avoid this fine tuning.) The gauge group is broken down to SU(3) 
X SU(2) X U(l) when S grows a VEV: (S) = Mdiag(2,2,2, -3, -3). The resulting super- 
heavy spectrum includes a heavy color Higgs chiral multiplet (3,1,|) of mass Mh = 5A2M, a 
vector multiplet (3,2,|) of mass Mv = 5V2gM, chiral multiplets (8,1,0) and (1,3,0) of mass 
Ms = |AiM, and a Standard Model gauge singlet chiral multiplet (1,1,0) of mass |AiM, 
where the numbers in parentheses are the SU(3) and SU(2) representations and hypercharge 
quantum numbers. 

In the following we assume 0.1 < A12 < 2.0, i.e. 10^'^ ^ Ci\i,2 ~ V^- The upper bound 
is imposed so that the model stays within the perturbative domain, while the lower bound 
excludes any anomalously small couplings. 

In addition to the renormalizable interactions, one may add the dominant non- 
renormalizable operator from Planck scale physics [|13|, 



where Mpi = 1/ =1/ VSvrG. In the first part of this paper we impose, for naturalness, 
|c| < 1. The main effect of this term is to modify the unification condition when S grows 
a VEV. Note for now that (S) will be O(Mgut) so this term enters with a coefficient 



Mqut/Mpi ~ (1/10 — 1/100) and we would naively expect its effects to be small JTl 



We will concentrate mainly on gauge coupling unification. The running of the gauge 
couplings with respect to the energy scale, /x, is governed by the two-loop renormalization 
group equations [|I5|, 

jt^) = - + (3) 

where ctj = gf /An, t = (In fi)/27T, fi = ^/ (arbitrary mass parameter), and hj are the Yukawa 
couplings. In MSSM, one-loop coefficients are hi = (33/5, 1, — 3)j. The two-loop coefficients 
are also well known and can be found elsewhere, e.g. in [0. 

We first discuss the effect of GUT scale thresholds without the Planck term. The GUT 
degrees of freedom are included in the running at their respective thresholds by 

ar\f^) = «.o-^(/i) - E A&L'^ln(^), (4) 

where ajo(/i)'s are calculated numerically to two-loop accuracy from their low energy values 
via the RGE's using the MSSM beta functions. The low energy values of aio(yu)'s include 
SUSY threshold which will be discussed later. In Eq. (^, the index a sums over the GUT 
degrees of freedom with masses less than fi, Abj^ = (0,2,3)j, A6*^ = (2/5,0,l)i, and 
A6*Y = (—10, —6, — 4)j. The largest of the Mq is called Mu, as this is where the coupling 
constants actually meet. Thus the unification condition is aj(Mu) = a5(Mu). 

On the other hand, if the Planck term is included, then when the VEV of S is inserted 
into the dominant Planck scale operator of Eq. (^, the kinetic terms for the gauge bosons 
will receive a contribution. Thus the unification condition will be modified by replacing 
a5(Mu)-^ by 

where M is the mass parameter entering in (S). 

At low energies we must consider the decoupling of supersymmetric degrees of freedom. 
This can be described at the one-loop level by three SUSY threshold parameters Mj, one 
for each coupling constant The meaning of these parameters is as follows: if above 



Mj we assume the threshold particles to be massless, and below Mj we assume them to 
be completely integrated out, and we assume the couplings meet smoothly at Mj, then at 



scales far from M, our running coupling constants will match the exact ones ||T^. Such an 
Mi can always be found for each i so long as the one-loop beta function above the threshold 
is different from that below |[T^. Thus the effect of SUSY thresholds is given by 



«r^(Mz) = a,o-^(Mz) + A6f ^ In(^). (6) 

Here, aj(Mz) are the couplings at Mz, while ajo(Mz) are the couplings one would obtain at 
Mz if one ran with the full SUSY beta function down to Mz. Abf = (5/2,25/6,4), gives 
the contribution to the f3 function from the additional SUSY degrees of freedom. This is 
sufficient for a two-loop analysis as well, so long as the SUSY thresholds are not too far 
from Mz. We treat all Standard Model degrees of freedom except the top as degenerate 
with or lighter than Mz. We take the top mass to be 174 GeV in accordance to the latest 
experimental indication . Uncertainties in its value do not affect our results significantly. 

There are two subtleties involved in relating the Mj's to the sparticle spectrum. First, 
the coupling constants do not actually jump in slope when the scale reaches the mass of a 
particle, rather they change gradually. If our measurements are performed in the region of 
changing slope, "match and run" will be inaccurate. Second, if the threshold region is close 
to the electroweak symmetry breaking scale, there are always mass splittings among the 
particles in the gauge multiplets. This potentially large effect has never been fully treated, 
and is essential for incorporating detailed SUSY spectra in unification analyses. Regardless 
of these subtleties, any supersymmetric spectrum can be accommodated by Eq. (|]) if the 
Mj's are allowed to vary below Mz as well as above [Q. 

Combining Eqs. (^), @, @, and the unification condition as modified by (|^), we arrive 
at the equation which we use for our calculations: 

1, , / M M M 

a5-^(Mu) 1-C-— ,l-3c-— ,l + 2c 



«-i(Mz)-6,— ln( 



Mpi Mpi Mpi 

1 - .Mu, 



27T 'Mz 



- iy.iM dt+T. f-./^^ jy,(i) 

In order to provide a bound relating to Mh and Mj we fix c and sin^(^vF) to take 
their maximum values [^| and Ai to take its minimum value, while A2 varies. At each 



point we iterate numerically in order to find the solution of Eq. (|^ and thus 03 and Mh- 
Some of the results are displayed in Fig. |I|, where we have chosen the degenerate case of 
Ml = M2 = M3 = MsusY, and we have plotted a bound each for MgusY = 10 GeV, 100 GeV, 
and 1000 GeV. 

The relevant bound on 03, the lower horizontal line in Fig. |l], can be parametrized by 



OtZ^rava = 0.040 + 0.0139 

-0.00579 tsusY - 0.00454 t^^gy, (8) 

where ta = (l/27r) ln(MQ,/Mz) for mass Mq,. For the case that the Mj's are non-degenerate, 
we may still use Eq. (||) to a close approximation, when we replace tsusY using the following 
formula: 

^susY = 0.305 ti + 7.738 ts - 7.043 ^3 
+2.38 + 20.91 + 6.90 tl 

-20.21 tit2 + 12.54 tit3 - 22.53 ta^s- (9) 

We find it useful to think of Msusy as defined in this formula as an effective mass scale to 
account for the SUSY thresholds in Eq. (P). Thus, in Fig. |l], although we have only plotted 
the results for the degenerate SUSY thresholds case, we may think of the bottom curves 
(which give the bound) as valid for all Mj with corresponding Msusy as given by Eq. 



Current bounds on the p ^ u + decay mode |Q imply Mh > 1.2 x 10^^ GeV 
Using Msusy = 100 GeV, (a characteristic value consistent with proton decay data), one 
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finds as > 0.112. Thus resolving the measurements is a crucial test of the model. 

We also plot in Fig. ^ the allowed region (shaded) where c = 0, and all GUT thresholds 
are taken to be degenerate (i.e. Mh = My = Ms = Mu), for the case of degenerate SUSY 
thresholds, and Mj between 10 GeV to 1000 GeV. We see that the lower bound on Mh 
implies now > 0.1145, and there is a strong correlation between a3(Mz) and Mh. 

We turn now to the Planck scale parameter c. Due to a cancellation, fixing sin^ {9w), 



a3(Mz), and the Mj's uniquely determines Mh when c = [^]. Assuming in the future we 
have a complete picture of the low energy physics, including the SUSY spectrum and the 
proton decay rate in both major modes, we can compare this value for c = (Mh(0)) with 
the correct Mh as determined from proton decay, and therefore determine c. The result can 
be summarized in the formula 

lOvrM Mh(0) ^ ^ 

M could be determined by the p + decay mode, since this decay is governed by 

My = 5^/8^Ta5 M. One may obtain from Eq. (|^). An approximate value is = 1/23. 
Note that a rough determination of c is possible even without the pion decay mode, e.g. 
by taking A2 = 1 so that M = Mh/5. In this approximation, the front factor in Eq. (|TUp 
is ~ 0.3, and while the ratio Mh/Mh(0) can vary substantially for reasonable range of low 
energy physics data, owing to the logarithmic dependence, c determined this way is in fact 
0(1) and not expected to be anomalously large. 

Although we have not yet observed proton decay, nor been able to determine the SUSY 
spectrum, using our naturalness conditions that 0.1 < Ai,2 < 2.0 and 10 GeV < Mi^2,3 < 
1000 GeV, and requiring Mh > 1.2 x 10^^ GeV we can already place bounds on 
c for different values of a3(Mz) and sin^^Ow)- We plot in Fig. ^ the allowed region for 
sin\ew) = 0.2294 and for sin^iOw) = 0.2327 in the c-a3(Mz) plane for Mi = M2 = M3. We 
find that in this analysis where Mj's are taken to be degenerate, the Planck scale term must 
exist for some values of siii^^Ow) and a3(Mz). We further note that ranging over the values 
of sin2(^^y) between 0.2294 and 0.2327, the LEP range for a3(Mz), i.e. 0.117 < < 0.129 



corresponds to —6.0 < c < 3.4; and the low energy data |^ range, 0.107 < 0:3 < 0.117, 
corresponds to —2.2 < c < 8.1. 

In contrast to the conclusions of Refs. IQ-Q , we have found that GUT models of the type 
considered here do have significant experimental consequences. Thus for models without 
Planck scale terms (c = 0), measurements made purely at the low energy (electroweak) scale 
can allow a prediction of the proton lifetime, and thus allow a direct test of the model. When 
c is left arbitrary, this is no longer possible. However, then low energy measurements will 
allow an experimental determination of c, and one has the remarkable possibility of seeing 
experimentally, for the first time, whether Planck scale physics exists. We note further, that 
the presence of GUT threshold effects and Planck scale terms do not qualitatively change 
the grand unification conclusions of Ref . . 

This work was supported in part by the National Science Foundation under grant number 
PHY-9411543. 
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FIGURES 

FIG. 1. Allowed regions in the parameter space projected onto the MH-a3(Mz) plane. The 
curves arc for MsusY = 10 GeV (solid), 100 GeV (dashed), and 1000 GeV (dot dashed). The 
shaded region corresponds to the case where c = 0, with all the GUT threshold mass scales taken 
to be degenerate, and MsusY allowed to vary between 10 GeV to 1000 GeV. 

FIG. 2. Allowed regions in the parameter space projected onto the as{Mz)-c plane. The curves 
are for sm'^{ew) = 0.2294 (solid) and sin'^{0w) = 0.2327 (dashed). 
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FIG. 1. Allowed regions in the parameter space projected onto the MH-a3(Mz) plane. The 
curves are for MsusY = 10 GeV (solid), 100 GeV (dashed), and 1000 GeV (dot dashed). The 
shaded region corresponds to the case where c = 0, with all the GUT threshold mass scales taken 
to be degenerate, and MsusY allowed to vary between 10 GeV to 1000 GeV. 
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